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Abstract 

A model of the Universe as a mixture of a scalar (inflaton or rolling tachyon from the string 
theory) and a matter field (classical particles) is analyzed. The particles are created at the 
expense of the gravitational energy through an irreversible process whereas the scalar field 
is supposed to interact only with itself and to be minimally coupled with the gravitational 
field. The irreversible processes of particle creation are related to the non-equilibrium pressure 
within the framework of the extended (causal or second-order) thermodynamic theory. The 
scalar field (inflaton or tachyon) is described by an exponential potential density added by a 
parameter which represents its asymptotic value and can be interpreted as the vacuum energy. 
This model can simulate three phases of the acceleration field of the Universe, namely, (a) 
an inflationary epoch with a positive acceleration followed by a decrease of the acceleration 
field towards zero, (b) a past decelerated period where the acceleration field decreases to a 
maximum negative value followed by an increase towards zero, and (c) a present accelerated 
epoch. For the energy densities there exist also three distinct epochs which begin with a scalar 
field dominated period followed by a matter field dominated epoch and coming back to a scalar 
field dominated phase. 

1 Introduction 

The solution of the problems of flatness, horizon and unwanted relics by the inflationary theory on 
the one hand and the recent measurements of the anisotropy of the cosmic microwave background 
and of the type la supernova SN 1997ff red-shift indicating that the Universe is flat with a present 
positive acceleration and a past decelerating period on the other hand, allows us to classify the 
evolution of the acceleration field of the Universe according to three major epochs, beginning with 
an accelerated inflationary period following a past decelerated epoch and leading back to a present 
accelerated phase. 

The early and present accelerated periods are dominated by a scalar field whereas the past 
decelerated phase is dominated by a matter field. In inflationary cosmology the scalar field is 
normally represented by the inflaton (see, e.g., the works QJI21EI) but recently a rolling tachyon 
- which has arisen from string theory |S] - was also considered as another candidate for the 
description of the accelerated phases of the Universe (see, e.g., the works El OD The 
inflationary period comes to an end when the scalar field has rolled to its potential minimum and 
begins to oscillate about it. By this time, the Universe is very cold and void. Some mechanism 
is then required to account for particle creation, that will fill and heat the Universe to allow for 
the standard hot Big-Bang evolution to take place. The most common mechanism adopted is the 
decay of the inflaton field in other particles, and it is known as reheating process. It goes without 
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1 



saying that this requires the coupling of the inflaton to other fields. However, the inflaton needs 
to be weakly coupled in order to inflation to work well. In view of that, in the present work 
we will make another approach and ignore the coupling of the scalar field to other fields. It is 
known from quantum field theory in curved space-time (see e.g. ^2)> that matter constituents 
may be produced quantum- mechanically in the framework of Einstein's equations. The energy of 
the produced particles is extracted from the gravitational field. Our procedure will then be to 
take this fact into account by a phenomenological description - first proposed, to the best of our 
knowledge, by Prigogine et al J3| ~ m which the particles are created during the evolution of the 
Universe at the expense of the gravitational energy through an irreversible process (see also the 
works [H[T£l[Tin[I7]). 

In the present work we have considered that the irreversible processes of particle production 
are related to a non-equilibrium pressure with an evolution equation coming from the extended 
(causal or second-order) thermodynamic theory (see e.g. ^3 El El El EH EDUCES) Moreover, 
since the scalar field potential does not require a minimum about which the inflaton will oscillate 
and decay in other particles, we shall adopt an exponential potential density added by a parameter 
which represents its asymptotic value and can be interpreted as the vacuum energy [231 1221 E3 HE! • 

We have shown, among other results, that a mixture of a scalar field (inflaton or tachyon) and a 
matter field can simulate the three phases of the Universe related to its acceleration field, namely, 
(a) an inflationary epoch with a positive acceleration followed by a decrease of the acceleration field 
towards zero, (b) a past decelerated period where the acceleration field decreases to a maximum 
negative value followed by an increase towards zero, and (c) a present accelerated epoch. For the 
energy densities there exist also three distinct epochs which begin with a scalar field dominated 
period followed by a matter field dominated epoch and coming back to a scalar field dominated 
phase. These results can also be obtained by choosing other types of potential densities for the 
scalar field (inflaton or tachyon). 

The work is organized as follows. In section 2 a system of three coupled differential equations 
for the scalar (inflaton or tachyon), acceleration and non-equilibrium pressure fields is determined. 
The solutions of the systems of coupled differential equations obtained in section 2 are found in 
section 3 for given initial conditions and for a given potential energy density of the scalar field. We 
close the work with a discussion of the results obtained in section 3. Units have been chosen so 
that c = h = k = 1. 



Let us consider a homogeneous, isotropic and spatially flat Universe modeled as a mixture of a 
scalar field and a matter field. The scalar field (inflaton or tachyon) represents a hypothetical 
particle while the matter field refers to the classical particles which are created at the expense of 
the gravitational energy. 

For this model of the Universe the energy-momentum tensor of the mixture is written as 



where the pressure and the energy density of the mixture are given in terms of the corresponding 
quantities for its constituents by p — p s +p m and p = p s + p m , with the indexes s and m denoting 
the scalar and the matter fields, respectively. We shall adopt the following convention: (a) s = <f) 
for the inflaton field and (b) s = ip for the tachyon field. Furthermore, (such that U^U^ — 1) is 

the four-velocity, denotes the metric tensor with signature (H ) whereas w refers to the 

non-equilibrium pressure. The non-equilibrium pressure is responsible for the irreversible processes 
of particle production [141 115} during the evolution of the Universe. 

The conservation law of the energy-momentum tensor ' ; „ = follows from Einstein's field 
equations and the Bianchi identities. In a comoving frame described by the Robertson- Walker 
metric it leads to the balance equation for the energy density of the mixture, i.e., 



2 Field Equations 



T flu = (p + p+ w)U^U u - + m)g' 



(1) 



p + 3H(p + p + w) = 0, 



(2) 
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where the quantity H = a(i)/a(t) denotes the Hubble parameter, a(t) is the cosmic scale factor 
and the over-dot refers to differentiation with respect to time t. 

The equation which connects the Hubble parameter with the energy density of the mixture is 
the Friedmann equation which - in a spatially flat Universe described by the Robertson- Walker 
metric - is written as 

H = — P, (3) 

where G is the gravitational constant. 
2.1 Inflaton field 

First we shall analyze the case where the scalar field is represented by the inflaton field </>(x^), 
which is described by the Lagrangian density 



C^-d^d^-V^), (4) 



where V{<j>) denotes the potential density of the inflaton field. The identification of the inflaton 
with a perfect fluid - i.e., by requiring that the inflaton interacts only with itself and it is minimally 
coupled with the gravitational field - allows us to write its energy-momentum tensor as 

= fa + p+WW - p+gT = fl'W ~ C*fT- (5) 

The last equality above is a consequence of Noether's theorem. 

By considering a homogeneous inflaton field and a comoving frame one can obtain from equations 
Q and (JSJ the relationships 



Ptt , = ±4> 2 + V(4>), 

V4 , = \4> 2 - Vtf), 



? - 2 (6) 



which connect the energy density p^ and the pressure p^, of the inflaton to its kinetic and potential 
energies. 

The time evolution equation of the inflaton field follows from the Euler-Lagrange equation which, 
in the homogeneous case, reads 

4> + 3H<j> + V'((f>) = 0. (7) 

Above, the prime denotes differentiation with respect to <f>. 

The evolution equation of the energy density of the inflaton field decouples from that of the 
matter field, since the differentiation of ©i with respect to time by taking into account the equation 
J7J leads to 

p 4 , + 3H(p 4> +p 4> ) = 0. (8) 
Hence the evolution equation of the energy density of the matter field can be written as 

p m + 3H(p m + p m ) = -3Hw. (9) 

thanks to J5J) and JSJ. We can interpret the term —3Hw in the above equation as the energy 
production rate of the matter field (see e.g., ^3 El E])- In a previous work ^H] one of the 
authors have calculated the energy-momentum pseudo-tensor of the gravitational field and found 
Xq° = —3H 2 /(8ttG). If we identify Xj£° with the energy density of the gravitational field pa, we 
can regard the term 3Hw as the energy production rate of the gravitational field: 

pc + 3H(pG - P4, - p m ) = 3i?tu, (10) 

that is, there is an irreversible energy flow from the gravitational field to matter creation. 

The matter field is supposed to obey a barotropic equation of state p m = w m p m where the 
coefficient w m may assume values in the range between < w m < 1. Some values for this coefficient 
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are: (a) w m = for dust or pressure- less fluid; (b) w m = 1/3 for radiation; (c) w m = 2/3 for non- 
relativistic matter and (d) w m — 1 for stiff matter or Zcl'dovich fluid. 

The equation which gives the time evolution of the cosmic scale factor can be obtained from 
differentiation of the Friedmann equation © with respect to time and elimination of p from the 
resulting equation, by using the balance equation for the energy density of the mixture (J2J). Hence, 
it follows 



H+-{w ri 



l]H 2 = AmG 



(w Tl 



i) y + K 



i)V(4>) 



(ii) 



thanks to ©, (j2Jl and to the barotropic equation of state of the matter field. 

Equation refers to a differential equation for the cosmic scale factor a(t) that depends on 
the non-equilibrium pressure zu(t) and on the potential density of the inflaton field V(<p). If we 
know a relationship between w{t) and a(t) for a given V((f>) it would be possible to find a solution of 
the system of differential equations Q and (| 1 1 1> for the time evolution of the cosmic scale factor and 
for the inflaton field. Here we shall consider that the non-equilibrium pressure obeys - within the 
framework of extended (causal or second-order) thermodynamic theory - the linearized evolution 



equation 



(12) 



Above, the coefficient of bulk viscosity r\ and the characteristic time r are considered as functions 
of the energy density of the mixture p, i.e., rj = ap, and r = rj/p where a is a constant (see 

e.g. H3GSEIU3). 

Now we have a system of three differential equations (JTJl , l|ll|) and (|12() , and in order to solve it 
we introduce the dimensionless quantities 



H = H/H 0l t = tH , w = m[8irG/(3H %)}, 

a = aH , V = V[8ttG/(3H§)], 4> = (/>y/8nG/3. 



(13) 



Above, the Hubble parameter Hq at t = (by adjusting clocks) is related to the energy density 



of the inflaton field at t — by H = y8irGp^/3, since we have assumed that at t = the energy 

density of the matter field vanishes (p°„ = 0) and 4>(0) = 0. 

With respect to the dimensionless quantities (|13|l the system of differential equations (JJJ, (flip 
and (|T2*)l reads 

'<j> + m<j) + V'{<l)) =0, (14) 



H+-(w n 



1)H Z = - 
' 2 



(w m - 1)— + (w n 



i)v(4>) 



zd + am = —3aH 



(15) 
(16) 

The time evolution of the inflaton field (f>(t), of the cosmic scale factor a(t) and of the non- 
equilibrium pressure vu(t) can be determined from the system of differential equations (I14fl through 
<|ltj|) once: (a) a potential density of the inflaton field V(<p) is chosen; (b) initial conditions at 
t = (by adjusting clocks) are given for the cosmic scale factor a(0) and its derivative a(0), for 
the inflaton field 0(0) and its derivative (f>(0) and for the non-equilibrium pressure vj(0); (c) values 
for coefficient w rn related to the barotropic equation of state of the matter field - and for the 
coefficient a - related to the irreversible processes of particle production - are selected. 

From the knowledge of the fields 4>(t), a(t) and w{t) it is possible to determine the time evolution 
of the energy densities of the inflaton and matter fields from 



1 ;r 



= H Z - 



1 v 



-y + v{4>) 



(17) 



thanks to ©i and (3J|. 



2 One is referred to 1271 for a derivation of the evolution equation 1121 within the framework of the Boltzmann 
equation. 



4 



2.2 Tachyon field 

Recently, Sen [5] has shown from a string theory that a rolling tachyon can be described by the 
Lagrangian density 

C v = -V{fp)y/\ - dptpdUp, (18) 

where V(<p) denotes the potential density of the tachyon field ip(x^). The corresponding Born-Infeld 
action for the tachyon field reads 



5' 



d X\J — g 



R 



167rG 



(19) 



In the above equation R represents the curvature scalar. 

From the action 11911 it follows the energy-momentum tensor of the tachyon field 

T$" = (p< P +p lp )U' i U v -p v g> a ', 
where the energy density p v , the pressure p v and the four- velocity are identified as 

V(<p) 

P v = -VQp)y/l - d^pd^p, 



(20) 



(21) 



y/ d v tpd"ip 

The field equation for the tachyon is obtained from the Euler-Lagrange equation and reads 



+ 



d»ipd v <p 

1 — daipd^if 



V'(<P) 
V(<p) 



0. 



(22) 



For a homogeneous tachyon field in a spatially flat Robertson- Walker metric, equations (|21[h 9 
and (|2"2"|l reduce to 



Pv 



V'(<p) 
3Hp + T7 ^-f = 0, 



(23) 
(24) 



1 - ip 2 r V(ip) 

respectively, whereas the four-velocity i|21|h becomes that of a comoving frame. 

The differentiation of equation IL'ol) 1 with respect to time leads to the evolution equation of the 
energy density of the tachyon field, i.e., 



p v + 3H(p v +p v ) = 0, 



(25) 



thanks to [12: il) and 1)24(1. Hence, the energy density of the matter field decouples from the corre- 
sponding equation for the tachyon field, has the same expression as that given by iJJJJ and can be 
interpreted in the same manner as in the previous section. 

Furthermore, from the differentiation of the Friedmann equation © with respect to time and by 
taking into account the barotropic equation of state for the matter field, it follows the dimensionless 
equation 



V(cp) 



= (w m + 1 - if ) 



(26) 



which relates the time evolution of the cosmic scale factor with the non-equilibrium pressure and 
the potential density of the tachyon field. The dimensionless quantities in this case are 



H = H/Ho, t = tH , zu = m[8-n:G/(3Hg)} 
a = aH , V=V[8ttG/(3H$)), <p = <pH 0l 



(27) 
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where the Hubble parameter Ho at t — (by adjusting clocks) is connected with the energy density 

p° v of the tachyon field at t = by H = ^JSnGp^/S. Here we have also assumed that at t = the 

energy density of the matter field vanishes (p^ = 0) and ip(0) = 0. 

The system of differential equations we have to solve now consists of the evolution equations: 
(a) for the non-equilibrium pressure l(Tp|l ; (b) for the tachyon field lt2"4"l) and (c) for the cosmic scale 
factor (|26[) . Once the time evolution of these fields is known one can determine the time evolution 
of the energy densities of the tachyon and matter fields from 

Pv V{<P) Pm „ 2 Vijf) / 08 % 
70" = Aj rf > If ^ ~ R 72' 

P v v 1 -^ P<p v 1 -^ 
3 Results and Discussions 

In order to solve the two systems of coupled differential equations for: 

(a) the non-interacting inflaton field [eqs. I|14|). I|15|l and iJTSJ], 

(b) the non- interacting tachyon field [eqs. (|16[) . (|24[1 and i|26[) ]. 

we have to choose at first the initial conditions at t = (by adjusting clocks). Here we have 
assumed: 

(i) a(0) = 1 for the dimensionless cosmic scale factor, 

(ii) H (0) = 1 for the dimensionless Hubble parameter, 

(iii) the inflaton field 0(0) and the tachyon field (p(0) were chosen in such a manner that the initial 
values of the corresponding potential densities were given by V(cj)(0)) = V(ip(Q)) = 1, 

(iv) (f)(0) = and (p(0) = so that the initial values of the energy densities of the inflaton and 
tachyon fields refer only to the corresponding potential densities, and 

(v) tu(0) = so that the irreversible processes of particle production begin just after the time 
t = 0. 

Among several models for the potential density of the scalar field (see, for example, Liddle and 
Lyth j28] ) we fix our attention to the exponential potential 

J V((f>) — exp[— p4>{t)\ + A, for the inflaton field, , , 

1 V(ip) = exp[— pip(t)] + A, for the tachyon field, 

where fi > and < A -^i 1 are free parameters, the first one is connected with the slope of the 
potential density and has influence on the velocity of the field which rolls down toward the potential 
minimum, whereas the second one represents the asymptotic value of potential density V and can 
be interpreted as the vacuum energy |231 1241 051 126| . 

Apart from the two free parameters A and y, there still remains much freedom to find the solutions 
of the four systems of coupled differential equations, since they do depend on the parameters: (a) 
w m which is related to the barotropic equation of state of the matter field, and (b) a which is 
connected to the transfer of energy from the gravitational field to the matter field. 

Before discussing the results in detail, let us briefly describe the evolution of the system. In the 
beginning all the energy is in the form of potential energy of the scalar field, which has negative 
pressure and therefore leads to an inflationary expansion of the Universe. As the scalar field begins 
to roll down its potential V , its kinetic energy increases and its potential energy decreases, raising 
its pressure and making its total energy decay. At the same time, the non-equilibrium pressure 
w increases in module and the irreversible processes of particle production begins. Therefore, the 
increasing production of matter and the decay of the potential energy of the scalar field contribute to 
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Figure 1: Acceleration a vs time t: inflaton field - straight line, tachyon field - dashed line. 



slow down the expansion of the Universe and the decelerated period begins. During this transition 
process the non-equilibrium pressure had started to decrease in module, and a little time after 
deceleration begins w was already driven towards zero. That means that the irreversible processes 
are only important during the early Universe, and there is no matter creation during the decelerated 
phase and thereafter, so that the matter fields evolve as usual. Moreover, the roll of the scalar field 
has been damped by the friction term (3H<p or 3H(p), such that its kinetic energy tends to zero 
and its potential energy tends to A. As long as p m 3> A the Universe evolves as in the standard hot 
Big-Bang model and the role of the scalar field is negligible. But when p m became comparable to 
A (the vacuum energy of the scalar field), a transition period takes place and the Universe begins 
to accelerate. 

The time evolution of the acceleration field a is plotted in figure 1 whereas in figure 2 are plotted 
the time evolution of the energy densities of inflaton p^, tachyon p v and matter p m fields. In these 
figures we have chosen the following values for the parameters: A = 0.005, p — 6.5, w m = 1/3, and 
a = 0.2. Below we shall comment how the changes of these parameters affect the solutions of the 
differential equations. 

We infer from figure 1 that there exist three distinct periods for the acceleration field, namely, (a) 
an inflationary epoch with an exponential growth to a maximum value followed by a decrease of the 
acceleration field towards zero, (b) a past decelerated period where the acceleration field decreases 
to a maximum negative value followed by an increase towards zero, and (c) a present accelerated 
epoch. We note from figure 1 that the acceleration field of the inflaton begins its decelerated period 
earlier than the corresponding one for the tachyon. This behavior is a consequence of the fact 
that the pressure of the inflaton field p^ = <f> 2 /2 — V{4>) may assume positive values, and in fact 
it does, due to the increase in the kinetic energy of the inflaton and the decrease in its potential 
energy. Hence, the inflaton field behaves as matter with positive pressure for a whil e, leadi ng to a 
precocious deceleration. By contrast, the pressure of the tachyon field p v = —V(tp)y/l — if 2 cannot 
assume positive values (for V(ip) > 0). Its contribution for the deceleration is not as intense as 
in the case of the inflaton field, since it can only come close to dust-like behavior (p v ~ 0). That 
prolongs the first accelerated period. 

We observe from figure 2 that for the energy densities there exist three distinct epochs which 
begins with a scalar field dominated period followed by a matter field dominated epoch and coming- 
back to a scalar field dominated phase. The following conclusions for the energy densities of 
the inflaton and tachyon fields can be obtained also from figure 2: (a) in the earliest times the 
energy densities of the inflaton and tachyon fields coincide, since <ft 2 <C V(<f>) and qp 2 <C 1 so that 
p^ pa V(4>) w V((p) w p v , (b) during the first accelerated phase the energy density of the tachyon 
field decays more slowly than that of the inflaton field, for two reasons: first, the influence of the 
downhill roll of the tachyon field (the increase in <p) on its respective energy density is stronger 
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Figure 2: Energy densities of infiaton p^, tachyon p v and matter p m fields vs time t: inflaton field 
- straight lines, tachyon field - dashed lines. 



than that of the inflaton field, since the former is proportional to 1/ \J\ — ip 2 , while the latter is 
proportional to <j> 2 ; second, the equation of motion of the tachyon field, (p = — (l — (p 2 )(3Hip+V'/V), 
tells us that the change in dp tends to zero as it approaches the value 1, therefore dp grows more than <j> 
and stands for a longer period nearby the value 1, and (c) in the transition from the past decelerated 
epoch to the present accelerated period the energy density of the tachyon becomes smaller than 
that of the inflaton, because the tachyon has rolled faster and further than the inflaton, implying 
that V(ip) < V{4>) and therefore p v < p^ - since in this period the potential energies dominate the 
kinetic energies. 

The behavior of the matter field can be interpreted as follows: the decrease in module of the 
non-equilibrium pressure vj for the tachyon field is slower - since its positive acceleration is more 
prolonged - and therefore more energy is transfered to the matter field than the corresponding case 
for the inflaton field. 

We shall comment upon the coefficient that is responsible for particle creation, namely a. By 
increasing the value of a: (a) the initial accelerated period grows due to the increase in module of 
the non-equilibrium pressure, and (b) the energy density of the matter field increases so that the 
matter dominated period becomes larger and the present accelerated period begins at later times. 

Let us now comment on the coefficient w m which refers to the barotropic equation of state of 
the matter field. By decreasing its value the positive pressure of the matter field decreases, and 
that implies: (a) a less pronounced deceleration, and (b) a slower decay of the energy density of 
the matter. 

If we decrease the value of the parameter A - which is related to the asymptotic value of the 
potential density V - we infer that: (a) there exists a longer matter dominated period since the 
vacuum energy A - which is the responsible for the present accelerated period - will dominate at 
later times, and (b) if A — > there exists no period of present acceleration, since p v w V(ip) — ► 
and p<j, w V{4>) — > 0. Hence, A is responsible for the present acceleration and can be interpreted as 
the dark energy. 

As was previously remarked the parameter p is related with the slope of the potential density 
of the scalar field. By decreasing the value of p the scalar fields <j> and ip roll more slowly so that 
the decay of their corresponding energy densities p<f, and p v is slower and the initial accelerated 
period becomes larger. Moreover, if p is very small the decelerated period for the inflaton may 
not exist due to the fact that the inflaton rolls more slowly when the slope of its potential density 
decreases. This behavior is not followed by the tachyon, since as was previously commented ip is 
pulled toward 1 and ip rolls more rapidly than <p. Hence, even for small values of p there exists a 
decelerated period for the tachyon field in this model. 
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As a final remark, we note that the same general behavior of the acceleration and of the energy 
densities fields may be obtained by choosing other usual types of potential densities for the scalar 
field found in the literature (see e.g. |28|1. 
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